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Multiple parton collisions will represent a rather common feature in pp collisions at the 
LHC, where regimes with very large momentum transfer may be studied and events rare in 
lower energy accelerators might occur with a significant rate. A reason of interest in large 
pt regimes is that, differently from low pt, evolution will induce correlations in x in the 
multiparton structure functions. We have estimated the cross section of multiple production 
of W bosons with equal sign, where the correlations in x induced by evolution are particularly 
relevant, and the cross section of bbhb production, where the effects of evolution are much 
smaller. Our result is that, in the case of multiple production of W bosons, the terms with 
correlations may represent a correction of the order of 40% of the cross sections, for pp 
collisions at 1 TeV cm. energy, and a correction of the order of 20% at 14 TeV. In the case 
of bb pairs the correction terms are of the order of 10 — 15% at 1 TeV and of the order of 5% 
at 14 TeV. 
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1. INTRODUCTION 

Multiple parton interactions in high energy hadronic collisions have been discussed long ago 
by several authors [1, 2]. Experimentally events with multiparton interactions have been first 
observed in pp collisions by the AFS Collaboration [3] and later, with sizably larger statistics, at 
Fermilab by the CDF Collaboration [4]. 

In multiple parton collisions the hadron is probed in different points contemporarily [2]. The 
non trivial feature of multiple parton collisions is hence its non-perturbative input, which has 
a direct relation with the correlations between partons in the hadron structure [5]. As the 
process is originated by the large population of partons in the initial state, the expectation 
is nevertheless that correlations should not represent a major feature in the process, with the 
exception of correlations in transverse space, which are directly measured by the cross section. 
Indeed the experimental analysis and most of the theoretical estimates have been done with this 
simplifying assumption and, although the statistics was too low to draw firm conclusions, the 
CDF Collaboration reported that the cross section is not influenced appreciably when changing 
the fractional momenta of initial state partons [4]. 

The much larger rates of multiple parton collisions expected at the LHC, with the possibility 
of testing different multiparton processes at different resolution scales, represents however a good 
motivation for reconsidering the approach to the problem. In particular an interesting process, 
where multiple parton collisions play an important role and which might be observed at the LHC, 
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is the production of multiple W bosons with equal sign [6], which would allow testing multiple 
parton interactions at a much larger resolution scale than usually considered. The evolution of 
the multiparton structure functions will play a non minor role in this case, leading to sizable 
correlations in fractional momenta. 

The purpose of the present note is to give some quantitative indication of the effects in 
multiparton collisions in a high resolution regime and compare with a case at a lower resolution. 
After recalling the basic features of the inclusive cross section of double parton collisions, we 
will hence evolve double parton distributions at high resolution scales. The effect of correlations 
induced by evolution will be estimated studying the cross sections of multiple production of equal 
sign W bosons and the cross section of multiple production of bb pairs, in the energy range 1 — 14 
TeV. 

2. DOUBLE PARTON CROSS SECTION 

With the only assumption of factorization of the two hard parton processes A and B, the 
inclusive cross section of a double parton-scattering process in a hadronic collision is expressed 
by [2, 7] 




where rij{xi, X2,b) are the double parton distribution functions, depending on the fractional 
momenta xi, X2 and on the relative transverse distance b of the two partons undergoing the hard 
processes A and B, the indices i and j refer to the different parton species and and a?i are 
the partonic cross sections. The dependence on the resolution scales is implicit in all quantities. 
The factor m/2 is a consequence of the symmetry of the expression for interchanging i and j; 
specifically m = 1 for indistinguishable parton processes and m = 2 for distinguishable parton 
processes. 

The double distributions rij{xi,X2,b) are the main reason of interest in multiparton collisions. 
The distributions (xi, 2:2, &) contain in fact all the information of probing the hadron in two 
different points contemporarily, though the hard processes A and B. 

The cross section for multiparton process is sizable when the flux of partons is large, namely 
at small x, and dies out quickly for larger values. Given the large parton flux one may hence 
expect that correlations in momentum fraction will not be a major effect and partons to be rather 
correlated in transverse space (as they must anyhow all belong to the same hadron). Neglecting 
the effect of parton correlations in x one writes 

Tij{xi,X2;b) = r,(xi)rj(x2)F](6), (2) 

where Tiix) are the usual one body parton distribution function and F'j{b) is a fimction nor- 
malized to one and representing the parton pair density in transverse space. The inclusive cross 
section hence simplifies to 

'^(A,B) = '^T.®kl^viA)a,i{B), (3) 

ijkl 

where aij{A) and dki{B) are the hadronic inclusive cross sections for the two partons labelled 
i and j to undergo the hard interaction labelled A and for two partons k and I to undergo the 
hard interaction labelled B\ 

0]^,= JcfbFlib)F^(b) (4) 
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are geometrical coefficients with dimension an inverse cross section and depending on the various 
parton processes. In the simphfied scheme above, the coefficients O^:'^ are the experimentally 
accessible quantities carrying the information of the parton correlations in transverse space. 

In the experimental search of multiple parton collisions the cross section has been further 
simplified assuming that the densities Fj do not depend on the indices i and j, which leads to 
the expression 

D m a{A)a{B) _ n 

where all information on the structure of the hadron in transverse space is summarized in the 
value of a single the scale factor, CTeff- In the experimental study of double parton collisions 
CDF quotes a^ff = 14.5m.b [4]. 

The experimental evidence is not inconsistent with the simplest hypothesis of neglecting cor- 
relations in momentum fractions, the resolution scale probed in the CDF experiment is however 
not very large, the transverse momenta of final state partons being of the order of 5 GeV. We 
will hence approach the problem in more general terms, focusing on multiple production of equal 
sign W bosons and of bbbb pairs, keeping into account the correlations in fractional momenta 
induced by evolution. 



3. TWO-BODY DISTRIBUTION FUNCTIONS 

The evolution of the double parton distribution function has been discussed in refs [8, 9] and 
more recently in [10]. The approach is essentially the same used to study particle correlations in 
the fragmentation functions [11], using the jet calculus rules [12]. 

Introducing the dimensionless variable 



27r6 



4 TT V M 



, 33 -2n^ 



where g'^i/J-^) is the running coupling constant at the reference scale /U^ and n/ the number of 
active flavors, the probability X2; t) to find two partons of types ji and j'2 with fractional 

momenta xi and X2 satisfy the generalized Lipatov-Altarelli-Parisi-Dokshitzer evolution equation 



where the subtraction terms are included in the evolution kernels P. 
If at the scale /j,^ one assumes the factorized form 

Di'^^zi,Z2,0) = Di\zi;0)Di'{z2;0)e{l -zi- Z2), (7) 

at a larger scale one obtains a solution which may be expressed as the sum of a factorized and 
of two non-factorized contributions: 

Di'^'{xi,X2;t) = Di'{xr,t) Di^X2;t)e{l-xi-X2) + Di]^^,{xi,X2;^ (8) 
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where the non-factorized contributions are expressed by the convolutions: 

''^ dzi dz2 



"&r,l<^l.»:2; *)=«(! -2^1 -12) 



■7-; J XI Jx2 ^2 Zi 



(9) 



,7:7-:, Jo ^xi ^1 Jx2 ^2 2l + ^2 



J J1J2 \zi+ Z2J ■'i ■'a Z2 



and the distribution functions Dj{x;t) satisfy the evolution equation 

dDj{x;t) 



dt 



Y.['^^A-':t)Pr-^{^). (11) 



with initial condition Dj{x]t = 0) = Sij 6(1 — x) 

Equations(ll) are solved by introducing the Mellin transforms 



Dj{n;t)= f\xx''Df{x;t), (12) 

which lead to a system of ordinary linear-differential equations at the first order. The solution 
is given by the inverse Mellin transform 

^^^^^ = I ^^^-''Di{n-t)=C-\Di{n-t),-Hx)), (13) 

where the integration runs along the imaginary axis at the right of all the n singularities, while 
represents the Inverse Laplace operator. 

The double distributions can then be obtained numerically. For inverting the Laplace Trans- 
form we have followed two different procedures [13]: the Gaver-Wynn-Rho (GWR) algorithm and 
the fixed Talbot (FT) method. The first procedure (GWR) is based on a special acceleration 
sequence of the Gaver functionals and requires to evaluate the transform only on the real axes; 
the second procedure (FT) is based on the deformation of the contour of the Bromwich inver- 
sion integral and requires complex arithmetic. Comparing the two methods we have found more 
stable results when using the (FT) method. The double distributions have hence been obtained 
by numerical integration with the Vegas algorithm [14], using the MRS99 [15] as input parton 
distribution function at the scale /i^. 

In the kinematical range of interest for the actual case (we never exceed x = .1) the contribu- 
tion of the term D^mrri ™ ^Q-(^) negligible. The first term in eq.(8) represents the factorized 
contribution usually considered and is the solution of the homogeneous (LAPD) evolution equa- 
tion, while the third term is a particular solution of the complete equation. 

The effect of the correlation terms induced by evolution is shown for gluon-gluon and for 
quark-quark in Fig. [1,2], where the ratio 

Rnn (X, X2-t) = ""^^rM^^-^^+^n^rM^--^^ .14. 

Di\x,,t)Di^{x2;t) 

is plotted as a function of x, with xi = X2 = x, with the following choice of parameters: 
/U = 1.2 Gey, rif = 4, factorization scale equal to the W mass, = 80.4 GeF (solid curves) 
and factorization scale equal to the bottom quark mass, m;, = 4.6 GeF (dashed curves). 
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As shown in Fig.[l], the ratio R^^ is nearly 35% for .t ~ 0.1 and decreases up to 8 — 10% for 
X ~ 0.01 and to 2% for x ~ 0.001, when the W mass is used as factorization scale. When taking 
the h quark mass as factorization scale, the value of the ratio is of the order of 10 — 12% for 
a; ~ 0.1 and decreases up to 5% and to 2% for x ~ 0.01 and x ~ 0.001 respectively. The ratio 
would of course be much larger (up to 60%) if going to larger x values. 

The ratio i?^"^ is shown in Fig. [2] for a few flavor choices. With the W mass as factorization 
scale, the ratios are of the order of 35, 20, 10% for x ~ 0.1, 0.01, 0.001. With the h quark mass 
as factorization scale the ratios are of order of 23, 10, 5% for x ~ 0.1, 0.01, 0.001. 
Apart from the case of hadron-nucleus collisions, when two different target nucleons take part to 
the process [16], the non-perturbative input of the double parton scattering cross section is not 
represented however by the distribution functions D-j^-'^{xi,X2;t) in Eq.(8), where all transverse 
variables have been integrated. The double parton scattering cross section, Eq.(l), depends in 
fact in a direct way also on the relative separation of partons in transverse space, which is of the 
order of the hadron size and hence outside the control of perturbation theory. 

Considering that the longitudinal and the transverse momenta of initial state partons are 
essentially decoupled in the process, because of the different scales involved, it's not unreasonable 
to assume phenomenologically a factorized dependence of the double distribution functions on the 
longitudinal and transverse degrees of freedom. Given the different origin of the terms in Dj^-'^ , 
it's also not unnatural to consider the possibility of having different non-perturbative scales, for 
the transverse separation of the factorized and of the correlated terms. In fact, although in the 
general case evolution would mix the two scales in the D-j^^^^^^ ^ term, the term D-j^^^^^^ ^ is very 
small in the kinematical regime of interest and the hypothesis of two different transverse scales 
is not inconsistent. 

We hence assume that the typical transverse distance between partons in -D^^j.^^^ and in 
^h cx>rr 1 corrcspouds to the relatively low resolution scale process observed by CDF and, to have 
an idea on the effects of the presence of two different scales in the double parton densities, we 
introduce a different transverse distance in the term D-j^^^^^.^. 2, related to the size of the gluon 
cloud of a valence quark, and corresponding to a relatively shorter range correlation term. The 
double parton distributions are hence expressed in the following way: 



While i^CTe// represents the transverse density of partons at a relatively low resolution scale, 
relevant in the kinematical conditions of the CDF experiment and leading to the measured value 
of the scale factor a^ff = 14.5 mb, iv^ is rather the transverse parton density characterizing 
partons correlated in fractional momentum, which becomes increasingly important when the 
resolution scale is large. To study the effect of the two scales we have let the smaller scale vary 
in the interval (7^0 '^eff assuming cTrg = 2.8 mb [17], which might represent the size of the 
gluon cloud of a valence quark in the hadron. To disentangle the effects of the correlation in 
fractional momenta we have neglected a possible dependence of the parton pair densities Fi(b) 
on the partons flavor. 



where the parton pair densities Fi{b) satisfy 




with i = O'eff} fr- 
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4. MULTIPLE PRODUCTION OF bb PAIRS AND OF EQUAL SIGN W BOSONS IN 

pp COLLISIONS 

For the purpose of the present analysis we have hence evaluated the contributions to multiple 
production of equal sign W bosons and to multiple production of bb pairs, due to multiple 
(disconnected) parton collision processes, taking into account the correlation terms in fractional 
momenta induced by evohition. 

As a matter of fact higher order corrections in as are very important in heavy quark produc- 
tion. To the purpose of the present analysis we have evaluated the cross section at the lowest 
order in perturbation theory, taking higher order corrections into account by rescaling the lowest 
order results with a constant factor K, defined as the ratio between the inclusive cross-section 
for bb production, cr(bb), and the result of the lowest-order calculation in pQCD. Our assumption 
is hence that higher order corrections in bbbb production may be taken into account by multi- 
plying the cross section of each connected process by the same factor K, so that higher order 
corrections are taken into account by multiplying the lowest order cross section by the iC-factor 
at the second power. In the actual calculation we have used a K factor equal to 5.7 ([18, 19]) 
and the value nib = 4:.6GeV for the mass of the bottom quark. 

The multiparton distributions have been obtained, as described in the previous paragraph, 
using as input distributions at the scale fi^ the MRS99 [15] parton distribution functions. Fac- 
torization and renormalization scale have been set equal to the transverse mass of the produced 
quarks. As for the dependence on the transverse variables, in addition to the usual factorized 
contribution, leading to the scale factor l/agff, in the present case the cross section includes 
also non factorized contributions, corresponding to the couplings of 2 both with -D^'j^ct 1 

and with 2- We have assumed a gaussian distribution for F^^^^{b) and for F„^{b). The 

scale factors are correspondingly 2/((Te// + cr^) and 1/(7^. 

In Fig.[l] we plot the gg correlation (the dominant contribution to 66 is gluon fusion) while the 
expected rise of the total 6666 cross-section is plotted in Fig. [3] {left-panel) as a function of the 
center of mass energy. The dashed curve refers to the double-parton scattering factorized term 
('^/act) given by eq.(5); the continuous curves refer to the double-parton scattering correlation 
contributions (o'^rr); with geometrical factors determined by setting r = tq (upper curve) and 
r = Teff (lower curve). The ratio between the contribution of the terms with correlations and 
the factorized term is shown in Fig. [3] {right panel) as a function of center of mass energy. The 
effect of the terms with correlations decreases by increasing the center of mass energy; depending 
on the values of r G [reff-ifo], correction effects may vary between (12 — 20)% at ^/s = ITeV and 
(3.5 — 6)% at \/s = lATeV . The decrease is faster as -^/i < 5TeV: for larger cm. energies the 
average fractional momentum < x > becomes smaller than 0.01, where the fraction R^^ stabilizes 
around 0.03 — 0.05, consistently with the amount of correction obtained for ^/s > 5TeV. 

In Fig. [4] we plot the 6666 production cross-section at -v/s = lATeV {left-panel) and at 
^/s = 5.5 TeV {right-panel), as a function of the minimum value of transverse momenta of 
the outgoing 6 quarks p™"^ in the pseudorapidity interval |r/| < 0.9. At ^/s = l^TeV with 
pmm g [0, 10]GeF one has < x >G [1.2,3.4] x lO"'"^, which leads to a contribution of the 
correlation terms of the order of (2 — 4)% and of (4 — 7)%, respectively for the lower and the 
higher choices of j?™" Fig.{o) (left-panel). At ^/s = 5.5TeV, in the considered range of variability 
Q^prnin Yyqs < X >€ [3.5, 6.5] X 10~^ and the contribution of the correlation terms can become 
of the order of 12%, ¥ \g.{5) {right-panel) . 

The cross sections of like-sign W pair production are evaluated at the leading order, hence 
including only quark initiated processes in the elementary interaction {qq' W). Higher order 
corrections are taken into account multiplying the lowest order cross section by the factor K ~ 



7 



1 + (87r/9) Q!s(^H/) [20]- We plot in fig. [6] (left-panel) the W^W^ cross-section as a function 
of the pp center of mass energy. As in the case of hbbb production, the dashed curve refers to 
the double-parton scattering factorized term {(y^act)i while the solid curves to the contribution 
of the terms with correlations {cTcorr)^ fo'^ the two different choices r = tq and r = re//. As 
one may infer from the behavior of the qq-correlation ratio, for < x >G [0.2, 6] X 10~^, which 
corresponds to the energy interval considered, the corrections due to the correlation terms range 
from (27 - 45)% at ^ = ITeV to (7.5 - 13)% at ^ = UTeV, depending on the choice of 
dr. The results for W~W~ production are presented in fig. (7). As shown in the right panel 
the correlation terms can give contributions ranging from (23 — 40)% at 1 TeV to (12 — 20)% at 
UTeV. 

5. CONCLUSIONS 

As an effect of evolution, the multiparton distributions functions are expected to become 
strongly correlated in momentum fraction at large and finite x [8-10]. On the other hand, 
the indications from the experimental observation of multiparton collisions at Fermilab [4] are 
not in favor of strong correlation effects in fractional momenta. The most likely reason being 
that the kinematical domain observed, relatively low x values and limited resolution scale, is far 
from the limiting case considered in QCD. 

The possibility of testing multiparton collisions at high resolution scales at the LHC will open 
the opportunity of testing the correlations predicted by evolution. To have an indication on the 
importance of the effects to be expected, we have considered a high resolution scale multiparton 
process (equal sign W pair production) and, for comparison, a sizably smaller resolution scale 
process {bbbb production) in pp collisions in the energy range ITeV < -^/i < 14Tev. In both cases 
the production process may take place either by single (connected) or by multiple (disconnected) 
hard parton collisions, while the two contributions may be disentangled applying proper cuts 
in the final state [4, 6, 19]. To study the effects of correlations we have hence worked out 
the disconnected contributions to the cross sections after evolving the multiparton distribution 
functions at high resolution scales. 

Our result is that the contribution of the terms with correlations, in equal sign W pairs 
production, might be almost 40%. of the cross section at 1 TeV and might still be a 20% effect 
at the LHC. The effect is much smaller in bbbb production, where corrections to the usually 
considered factorized distribution are typically between 5 and 10%. 
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Figure 2: qq correlation ratios for xi = X2 = x, with factorization scale equal to the W mass (sohd- 
curves) and to the b mass (dashed-curves). 
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Figure 3: Left-panel: bbbb inclusive cross-section as a function of center of mass energy; the dashed 
curve refers to the double-parton scattering factorizcd term cr^acf while the upper and the lower solid 
curves refer to contributions with correlations cr^^^, with geometrical factors obtained setting r = tq and 
r = reff respectively. 

Right-panel: ratio between cr^^^ and f^^^ as a function of center of mass energy. 
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Figure 4: bbbb production cross-section at \/i = lATeV and ^/s = 5.5 TeV as a function of p™*", the 
minimum value of transverse momenta of the outgoing b quarks; the dashed curves refer to the double- 
parton scattering factorized term a^^^^ , while upper and lower solid curves refer to (T^^r ^ith geometrical 
factors obtained setting r = ro and r = Veff respectively. 
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Figure 5: ratio between <t^^^ and cr^^^ at \/i = lATeV [left-panel) and \/i = 5.5TeV (right-panel) as 
a function of p™'". 
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Figure 6: Left-panel: W'^W^ inclusive cross-section as a function of the center of mass energy; the 
dashed curve refers to the double-parton scattering factorized term Uj^^j, while upper and lower sohd 
curves refer to u^^^, with geometrical factors obtained setting r = ro and r = Veff respectively. 
Right-panel: ratio between cr^^^ and cr^act ^ function of center of mass energy. 
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Figure 7: Left-panel: W^W~ inclusive cross-section as a function of the center of mass energy; the 

1 term crj^j^^, while the upper and the lower 
, with geometrical factors obtained setting 



dashed curve refers to the double-parton scattering factorizcd term crj^j^^ , while the upper and the lower 



solid curves refer to the contribution with correlations a^^^ 
r = ro and r = Veff respectively. 

Right-panel: ratio between cr^j,^ and crfact ^ function of center of mass energy. 
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